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ABSTRACT 

Context. The dominant astrophysical production site of the r-process elements has not yet been unambiguously identified. The suggested main 
r-process sites are core-collapse supernovae and merging neutron stars. 

Aims. We explore the problem of the production site of Eu. We also use the information present in the observed spread in the Eu abundances 
in the early Galaxy, and not only its average trend. Moreover, we extend our investigations to other heavy elements (Ba, Sr, Rb, Zr) to provide 
additional constraints on our results. 

Methods. We adopt a stochastic chemical evolution model that takes inhomogeneous mixing into account. The adopted yields of Eu from 
merging neutron stars and from core-collapse supernovae are those that are able to explain the average [Eu/Fe]-[Fe/H] trend observed for solar 
neighbourhood stars, the solar abundance of Eu, and the present-day abundance gradient of Eu along the Galactic disc in the framework of a 
well-tested homogeneous model for the chemical evolution of the Milky Way. Rb, Sr, Zr, and Ba are produced by both the s- and r-processes. 
The r-process yields were obtained by scaling the Eu yields described above according to the abundance ratios observed in r-process rich stars. 
The s-process contribution by spinstars is the same as in our previous papers. 

Results. Neutron star binaries that merge in less than 10 Myr or neutron star mergers combined with a source of r-process generated by massive 
stars can explain the spread of [Eu/Fe] in the Galactic halo. The combination of r-process production by neutron star mergers and s-process 
production by spinstars is able to reproduce the available observational data for Sr, Zr, and Ba. We also show the first predictions for Rb in the 
Galactic halo. 

Conclusions. We confirm previous results that either neutron star mergers on a very short timescale or both neutron star mergers and at least 
a fraction of Type II supernovae have contributed to the synthesis of Eu in the Galaxy. The r-process production of Sr, Zr, and Ba by neutron 
star mergers - complemented by an s-process production by spinstars - provide results that are compatible with our previous findings based on 
other r-process sites. We critically discuss the weak and strong points of both neutron star merging and supernova scenarios for producing Eu 
and eventually suggest that the best solution is probably a mixed one in which both sources produce Eu. In fact, this scenario reproduces the 
scatter observed in all the studied elements better. 

Key words. Galaxy: evolution - Galaxy: halo - stars: abundances - stars: neutron - stars: rotation - nuclear reactions, nucleosynthesis, 
abundances 


1. Introduction 

The heavy element Eu is an r-process element; that is to say, 
it is produced by neutron captures on heavy elements (for in¬ 


stance as Ee) in a rapid process, where rapid refers to the 
timescale of the neutron capture rates relative to the yS-decay 
rates of unstable n uclei. The main producti on site of Eu is still 
a matter of debate ( Thielemann et ahllzOl ll e. g.) 
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Observations of heavy element abundances in Galactic halo 
stars provide important constraints on the astrophysical site(s) 
of r-process nucleosynthesis. Interestingly, a wide spread is 
found in the [Eu/Fe] ratios in halo stars (as well as for sev¬ 
eral s-process elements, such as Ba, Y, La), and is much wider 
than the spread found for [tr/Fe] ratios in the same stars. 

Using an homogene o us ch emical evolution model, 
Mennekens & VanbeverenI ( 20141) conclude that neutron 


star/black hole mergers could be responsible for the Galactic 
r-process production. Their model did not consider a further 
contribution fro n t Type II supernovae (SNell). More recently, 
Matteucci et all (l2Q14l) have employe d a detailed chemical 
evolution model dRomano et al.l lioioh to study the evolution 
of Eu in the Galaxy. Two possibilities for Eu production were 
considered: i) production by core-collapse SNe (stars with 
initial masses from 9 to 50 M©) during the explosion, and ii) 
neutron star mergers (NSM). The classical production site for 
r-process elements, henc e Eu, is core-collapse SNe dTruranl 
198ll:ICowan et al.|[l99lh . 

The reason for introducing neutron star mergers as an 
alternative to Eu production resides in the large uncertain¬ 
ties present in hydrodynamical nucleosynthesis calculations 
for r-process elements in massive stars, in particular, that 
neutrino winds in SNII explosions are proton-rich or only 
slightly neutron-rich (see for example Arcones et al.l 12007 : 
Wanajo et all 1201 ll: lArcones & Thielemann 20131 and refer¬ 
ences therein) and therefore have difficulty producing Eu. 
At the present time, only the magneto-rotational driven SNe 
(MRD SNe) scena rio has been shown to be a promising source 
of the r-process bv lWinteler et al.l d2012l) in the context of mas- 
sive stars, and this resu lt has been conhrmed very recently by 
Nishimura et all d2015 ). However, given the specihc conhgu- 
ration needed by these progenitors, they are expected to be 
rare. On the other hand, the nucleosynthesis calculations rel¬ 
ative to neutron star mergers have provided robust results con¬ 
cerning the r-process element production in these objects (see 
for example Rosswog et al. 1999 . 2000 ; Oechslin et al. 2007 ; 


Bauswein et all 2013[ Rosswog 20131: iHotokezaka et al. 2013 ; 

Kvutoku et al.ll2013 ). 

It has been suggested that up to 10“^ M© of r-process matter 
may b e ejected in a single coalescence event. iMatteucci et al.l 
d2014l) suggest that neutron star mergers could be entirely re¬ 
sponsible for the Eu production in the Galaxy if the coalescence 
timescale is no longer than 1 Myr for the bulk of neutron star 
binary systems, the average Eu yield is 5 • 10“^ M©, and the 
mass range of progenitors of neutron stars is 9-50 M©. They 
also conclude that a mixed scenario could be acceptable, where 
both merging neutron stars and core-collapse SNe contribute to 
the Eu production. In the mixed scenario, the Eu yields from 
merging neutron stars should be lower since core-collapse SNe 
contribute to the enrichment. In particular, it was concluded 
that SNe in the range 20-50 M© should produce 10“^-10“^ M© 
of Eu each. Both models could reproduce the average trend of 
[Eu/Fe] versus [Fe/H] in the solar neighbourhood, the solar Eu 
abundance, and the Eu abundance gradient along the Galactic 
disk. 

By relaxing the instantaneous mixing approximation, it is 
also possible to explore the information contained in the ob¬ 


served scatter (or lack of ) in the different abu ndance ratios and, 
in particular, in [Eu/Fe]. lArgast et al.l (120041) explored the im¬ 
pact of Eu production by merging neutron stars and SNell In 
their model the diffusion of the stellar ejecta into the interstellar 
medium is treated dynamically, hence predicting the chemical 
spread in the c hemistry of the surroun d ing gas. Unlike the ap¬ 
proach taken in IMatteucci et al.l (l2014l) . lArgast et al.l (l2004 do 
not provide a model where both the NSM and core-collapse 
SN are simultan eously taken in to account. Another approach 
was presented bv ICescuttil (l2008h who modelled the inhomoge¬ 
neous mixing of the interstellar medium by means of a stochas¬ 
tic chemical evolution model. However, in this case the au¬ 
thor considered only SNell a s a site of production of the r- 
process. Both in lArgast et^ (l2004l) and ICescuttil (l2008h . a 
large abundance scatter is predicted for neutron capture ele¬ 
ments, whereas a much lower scatter is found for alpha/Fe 
abundance ratios, in agreement with observations. In both 
cases this was interpreted as a consequence of the stochastic 
formation of massive stars coupled with the different stellar 
mass ranges fro m which different elements come. In particular, 
Cescuttil ( 2008h suggested that the wide spread observed in neu¬ 


tron capture elements and the signihcantly narrower spread in 
Qf-elements occurs because the site of production of a-elements 
includes the whole range of massive stars from 10 to 80 M© 
whereas the mass range of production for neutron capture ele¬ 
ments lies between 12 and 30 M©. More recently, cosmologi¬ 
cal SPH simulations that include treating the chemical elements 
have investigated the NSMs as possible sou rces of the r-process 
dvan de Voort et al.ll2015l:IShen et al.ll2014l) . 

It is now impor tant to check for consi stency between the 
results obtained by IMatteucci et al.l (l2014l) on Eu with other 
r-process elements as well. In particular, it has recently been 
shown that magneto -rotationahy driven (MRD) supernovae 
( Winteler et al .11201 ^ represent a promising sourc e of r-process 
in the early Galaxy dCescutti & Chiappinill2014l) . This model 
was able to reproduce the observed spread in the abundance ra¬ 
tios not only of Eu, but also of Sr, Ba, and Y. In the case of 
Sr, Ba an d Y parts of the production most likely came from 
spinstars dPignatari et al.|[200^: iFrischknecht et ^12012 ). and 
the spinstar contribution to Eu is expected to be negligible. 
It has been shown that including the contribution of spinstars 
plays a key role in explaining the long-standing problem of the 
observed scatter in [Sr/Ba] in t he Galactic halo, as hrst sug¬ 
gested in IChiappini et al.l d2Qllh and lat er demonstrate d by th e 
inhomogeneous model calculations of ICescutti et al.l d2Q13 ). 
Interestingly, the spinstar scenario also plays a key role in 
explaining light element observations such as C and N (see 


IChiappini et aLlbOO^ l2008l: ICescutti & Chiappinill201C)l) . 

The goal of the present work is to evaluate the impact of 
including a neutron star merging scenario (which produces Eu, 
but also the other n-capture elements mentioned above) on the 
previous conclusions based on the MRD SNe plus spinstar sce¬ 
nario. In the present work we again provide our predictions 
for Ba and Sr but also include two other light neutron-capture 
elements that have not been modelled before, Zr and Rb. In 
parti cular, Rb has only been measured in a few globular clus¬ 


ters (Barbuv et al.ll20l4 Yong et al.ll20l4lD’Orazi et al.|[2Q13 

Yong et al.ll2008;lwaherstein et al.ll20Q7h. 
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The paper is organized as follows. In Section 2 we summa¬ 
rize the observational data considered in this work. In Section 
3, we introduce the chemical evolution model. In Section 4 we 
present our results and compare them to the available observa¬ 
tions. In Section 5, we draw some conclusions. 

2. Observational data 


dances of their heav y elements show fe atures compatible to 
those of normal stars ( Cescutti et al.ll2013 ). Nevertheless, since 
they are not the main focus of this work, CEMP-no stars are not 
included in the present analysis. We did not show any data for 
Rb, since for this element we could not find data for field stars 
at extremely low metallicity, and to our knowledge the lowest 
metallicity data for Rb in the literature refer to stars with [Fe/H] 


http://cdsarc.u-strasbg.fr/cgi-bin/qcat7J/AN/331/474 


^ The list of authors we use from the collection are McWilliam et al. 

(1995 

). iMcWilliaml (Il998l). 

Westin et al.1 (l2000h. lAoki et al. 

(2002 

). ICowan et al.l (I 2 OO 2 I). 

Ivans et al.l (l2003l). iHonda et al. 

(2004 

). lAoki et al.l (l2005l). Barklem et al.l (l2005l). lAoki et al. 

(2006 

). Ilvans et al.1 (l2006h. Masseron et al.l (I 2 OO 6 I). Preston et al. 

(2006 

). Aoki et al.l (l2007b. iFrancois et al.l (l2007b. iLai et al.l ( 

2007), 

Cohen et al.l ( 

20081). iLai et al.l ( 

20081). iRoederer et al.l ( 

2008), 

Bonifacio et al. 

(l2009h.lHavek et al. 

(|2009|) 


We employed the sa me data as us ed in ICescutti et al.l (120131) : 
the data compiled bv iFrebell (l201Qli1 and labelled as halo star£|. 
We excluded all upper limits and carbon-enhanced, metal-poor 
(C EMP) stars. For CEM P stars we adopt the definition given 
bv iMasseron et aP (I 2 OIQI) . where a CEMP star has [C/Fe]>0.9. 
In the compilation by Frebel, there is also a large portion of 
stars without carbon measurements. For these stars we can¬ 
not establish whether they are CEMP stars or not; neverthe¬ 
less, since they represent a large portion, we decided to include 
them in our plots, but to distinguish them graphically from 
the confirmed normal stars. CEMP-s stars are excluded from 
our comparison because, if carbon enhancement is caused by 
transfer of matter from an evolved companion, the abundances 
of s-process elements are likely to be affected, too. Therefore, 
no meaningful comparison can be done with the predictions of 
our chemical evolution model in this case, since it refers to the 
chemical composition of the stars at birth. 

We expect CEMP-no stars to behave differently from nor- 
mal stars only in the abundances of their lig ht elements 
( Cescutti & Chiappinil I 2 OIQI: iMaeder et all 12014) . The abun¬ 


the dimensions of the volume are expected to be large enough 
to allow us to neglect the interactions between different vol¬ 
umes, at least as a first approximation. For typical ISM densi¬ 
ties, a supernova remnant becomes indistinguishable from the 
ISM - that is, merges with the ISM - before reaching ~ 50 pc 
( Thornton et al.|[l9^ : therefore, we decided to have a typical 
volume with a radius of roughly 90 pc, and the number of as¬ 
sumed volumes is 100 to ensure good statistical results. We did 
not use larger volumes because we would lose the stochastic- 
ity we are looking for; in fact, larger volumes produce more 
homogeneous results. 

In each region, we assumed the same law for the infall of 
the gas with pr imordial compo s ition, following the homoge¬ 
neous model bv IChiappini et al.l ( 2008h : 


dGaSinit) 
dt 


oc e 


-(t-toficrl 


( 1 ) 


where to is set to 100 Myr, and is 50 Myr. Similarly, the star 
formation rate (SFR) is defined as 


SFR(t) oc (pgasiO) 


1.5 


( 2 ) 


where Pgas(t) is the density of the gas inside the volume under 
consideration. Moreover, the model takes an outflow from the 
system into account: 


dGaSputit) 

dt 


ocSFR(t). 


( 3 ) 


-2 that be long to globular clusters measured by IVong et al 

(I2nn6ll2008h . 


3. The chemical evolution model 

The chemical evolut ion r nodel adopted here is the sa me as in 
Cescutti et all (l2013h and ICescutti & Chiappinil (l2014l) . We re¬ 
view its main characteristics here for the reader’s convenience. 

We considered the chem ical evolution model presented in 
Cescutti & Chiappinil ( 2Q1(^. which is base d on the inhomoge- 
nous model devel oped bv ICescu^ ( 2008b and on the homo¬ 
geneous model of IChiappini et al.l (l2008h . The halo consists of 
many independent regions, each with the same typical volume, 
and each region does not interact with the others. Accordingly, 


Knowing the mass that is transformed into stars in a time 
step (hereafter, we assigned the mass to one star with 

a random functio n, weighted according to the initial mass 
function (IMF) of IScalol J 19861) in the range between 0.1 and 
100 Mq. We then extracted the mass of another star and re¬ 
peated this cycle until the total mass of newly formed stars ex¬ 
ceeded In this way, is the same in each region 

at each time step, but the total number and mass distribution 
of the stars are different. We thus know the mass of each star 
contained in each region, when it is born and when it will die, 
assuming the stellar lifetimes of iMaeder & Mevn^ ( 1989h . At 
the end of its lifetime, each star enriches the ISM with its newly 
produced chemical elements and with the elements locked in 
that star when it was formed, excluding the fractions of the el¬ 
ements that are permanently locked in to the remnant. 

As shown in ICescutti et al.l (l2013h . our m odel is ab l e to re - 
produce the MDF measured for the halo bv iLi et akl ( 2Qld) . 
This comparison shows that the timescale of enrichment of the 
model is compatible with that of the halo stars in the solar 
vicinity. Moreover, our model predicts a small spread for the 
Qf-elements Ca and Si, which is compatible with the observa¬ 
tional data. 

3.1. Stellar yields for Eu 

For the Eu production sites, we consider NSMs and core¬ 
collapse SNe, as mentioned in the Introduction. To take the Eu 
production from NS Ms into account, we nee d to define the fol¬ 
lowing quantities (^see lMatteucci et al1l2014 ): 
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1. the fraction of massive stars belonging to double NS sys¬ 
tems that will eventually merge, namely the realization 
probability for such events, a^s', 

2. the time delay between the formation of the double neutron 
star system and the merging event, A%s; 

3. the amount of Eu produced during the merging event, . 


Concerning NSM yiel ds, we also follow what is assumed in 


Matteucci et al 



Korobkin et al 
from 10“^ to 10“^ M( 


= m;^“(0.01 + 1.98 • Rand(n)), 


jEu 


A^ns, lArgast et alJ (120041) and iMatteucci et alJ (120141) consid 


Cowan etal. 1991: Wooslev et al. 1994 

shimaru & Wanaio 

19991: Travaelio et a 

. I 999 I: Wanaio et al. 

200 it Arsastetal. 

2004 ICescutti et all 

20061). The yields of Eu from SNe II that 


production for the 50 M© to balance the higher production 
of iron in our models. Since the most recent results concern¬ 
ing the production of Eu in SNe do not confirm these chan¬ 
nels (see Arcones et al. 2007, for the high mass channel, and 
Wanajo et al. 2011, for the low mass channel), we investi¬ 
gate another possible production connected to massive stars: 


the ma gneto-rotational driven (MRD) scenario. IWinteler et al 


; in particular, we assume the yields of 
, who suggest that NSM can produce 
0 of Eu per event. 


In the present paper, we assume that a fixed fraction of 
all the massive stars formed in our simulation is a progenitor 
of NSMs and produces r-process material. The progenitors are 
randomly chosen among the massive stars formed in the range 
9-30 M©. The prescriptions for the different models are sum¬ 
marized in Table [U where we list (i) the model name (column 
1), (ii) the delay time for coalescence of NS in binary systems 
(column 2), (iii) the fraction of massive stars that are hosted 
in binary systems leading to NSMs (column 3); (iv) the mass 
ejected as newly produced Eu in NSM events (column 4), and 
(v) the Eu yields from massive stars. Eor Models NS03 and 
NS04, the amount of r-process in the single event is not con¬ 
stant, but we considered the possibility that the amount of mass 
ejected as r-process varies. Since the variation is unknown, we 
assumed the following range: the minimum production is 1% 
of the average Eu amount, and the maximum is twice the aver¬ 
age. Since the total production should be conserved, the ejected 
mass for the ^-th star (r-process producer) in these models can 
be described by the following equation: 


(4) 


where Rand{n) is a uniform random distribution in the range 
[0,11. This function is the same as the one presented in 


Cescutti & Chiappinil (120141) . In Table \T\ we report the aver¬ 
age Eu amount for the merging event (M^^) assumed for each 
model. 

E or the time delay du e to t he coalescence of the t wo NSs, 


ered different timescales: 1 Myr, 10 Myr, and 100 Myr. Here 
we consider the same timescales. It is worth noting that in pre¬ 
vious works, as well as in this one, it is assumed that all neutron 
star binaries have the same coalescence timescale. Clearly, a 
more realistic approach should consider a distribution function 
of such timescales, in analogy with SNela for which a dis tribu- 
tion for the explosion times is defined ( see lGreg^l2005 ). 


(12012[) have shown that the combination of high rotation and 
strong magnetic field in the inner core of a n exploding SN 
prom otes an r-process production (see also iNishimura et 
2015h . This specific configuration is rare, therefore only a 
limited number of S Ne have this fate, and as mentioned in 
Winteler et al. (12012h. it shoul d be more frequent in the metal- 
poor regime dVoon et al.ll2006h . Therefore, we explore a model 
in which 10% of the SNe produ ce r-process material, with th e 
same prescriptions as assumed in ICescutti & Chiappinil (12014 ). 
However, unlike that paper, we consider here that this channel 
is only active at low metallicity, Z < 10“^. Coupled with the 
production of MRD SNe, we also assume a NSM production 
with a fixed delay of 100 Myr (see Table [T]). 

Einally, we note that in the spinstar framework, Eu is pro¬ 
duced in negligible amounts. 

3.2. Stellar yields for Rb, Sr, Zr, and Ba 

Rb, Sr, Zr, and Ba are produced by both the s- and r- 
processes. The r-process yields are obtained by scaling the Eu 
yields adopted here acc ording to the abunda nce ratios observed 
in r-process-rich stars dSneden et aklbOOSl) . Another possible 
choice would be to take the solar system r-proc ess contribution 
into y cou nt as determined, for e xample, by ISimmerer et al.l 
(I2OO4I) and I Arlandini et al.l ( 1999h : we checked that this does 
not seriously affect the results for these elements. Our choice 
does not rest on the results of the theoretical computations 
of the main r-process; rather, we infer the r-process elements 
yields from an observational signature that could reflect a com¬ 
bination of processes (e. g. main r-process -r weak r-process). 
Moreover, a given process might present some intrinsic varia¬ 
tion. Indeed, from an observational point of view, it seems that 
the robust pattern for the r-process does not extend to the ele¬ 
ments between the first and second r-process pe aks even within 
the c lass of highly r-process-enhanced stars dRoederer et al 


Among core-collapse SNe, different candidates for Eu pro¬ 
duction have been studied in the past, and they can either 


we adopt here are similar to thos e of lArgast et alJ (1200^ (their 
model SN2050) modified as in IMatteucci et al.l ( 20141) (their 
model Mod2SNNS), as shown in Table [U They are coupled 
with a production from NSM with a delay of 1 Myr (see Table 
[I])- Actually, sin c e we a dopt iron yields from SNell larger than 
iMatteucci et all d2Q14l) . we need to slightly increase the Eu 


20141) . In the future, we plan to use theoretical results for the 
r-process ratio to investigate this aspect. 

The spinstars’ s-process c ontribution for all ou r models is 
the same as in the/^'-model of ICescutti et al.l d2013l) . However, 
we show here results for rubidium and zirconium, which were 
not treated in the previous paper, therefore we recall that for 
the yields at Z =10~^, we conside red the stellar yields ob¬ 
tained by lErischknecht et al.l ( 20121) after decreasing the re¬ 
action rate for ^^0(a,y) from ICaughlan & Eowleil (1988) by 
a factor of 10, which enhances the s-process productiorlj. 
Unfortunately, we only have results with this reaction rate for 
a single mass (25 M©) at Z=10“^, and we used the scaling fac¬ 
tor obtained for the whole range of masses (for more details. 


^ A value of the metallicity of Z= 10 ^ corresponds to [Fe/H] ^ -3.5, 
with small variations due to the stochasticity of the models. 
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Table 1. Prescriptions for Eu production for different models. See text for explanations. 


Model name 

ANs 


neutron star mergers 

^newlyproduced ) 

massive stars 

NSOO 

1 Myr 

0.02 

5-10"^ (constant per merging event) 

no production 

NSOl 

10 Myr 

0.02 

5-10"^ (constant per merging event) 

no production 

NS02 

100 Myr 

0.02 

5-10"^ (constant per merging event) 

no production 

NS03 

1 Myr 

0.02 

on average 5-10“^ (varying as in equation |4]) 

no production 

NS04 

1 Myr 

0.04 

on average 2.5-10“^ (varying as in equation |4]) 

no production 

NS+SN 

1 Myr 

0.02 

3-10"^ (constant per merging event) 

2- lO-’^-5-lO-’ 

(linear interp. in the range 20-50 M©) 

NS+MRD 

100 Myr 

0.02 

1.5-10“^ (constant per merging event) 

on average MO"^ for 10% of stars in S-SOM© 
(varying as in equation |4]) 


see 


Cescutti et al .1120131) . Indeed, there are no nucleosynthesis 


calculations for spinstars currently carried out with a reduced 
value of the y) rate for a metallicity higher than Z=10“^, 

an d we adopted those comput ed with the standard value given 
bv ICaughlan & Fowled (Il988b . We need to keep this caveat in 
mind when interpreting our theoretical predictions for the in¬ 
termediate metallicity range. We also considered the s-process 
contribution from sta rs in the mass range I.S-SM q, by imple¬ 
menting the yields bv ICristallo et al.l (l2009[l2Qlll) in the mod¬ 
els. We underline, however, that this production channel affects 
the model results only at moderate metallicity ([Fe/H] —1.5). 


4. Results 

4.1. NSM models 

In Fig [T] (left panel) we show the distribution of synthetic halo 
stars in the [Eu/Fe]-[Fe/H] plane as predicted by our inhomo¬ 
geneous model when assuming that (i) Eu is produced only in 
NSMs with progenitors for neutron stars in the range 9-30Mo,; 
(ii) 2% of massive stars are in binary systems with the right 
characteristics to lead to merging NS, (iii) each merging event 
produces 5-10“^ M© of Eu, and (iv) there is a fixed delay time 
of 1 Myr between the formation of the double NS system and 
the merging event (model NSOO, see Table [T]). 

It is seen that, while explaining very well the data of stars 
with [Fe/H]> -2.2 dex, such a model fails to explain the pres¬ 
ence of stars with [Eu/Fe]<0.5 for [Fe/H]< -2.5 and does 
not explain the existence of any star with Eu measurements at 
metallicities lower than [Fe/H]= -3. The upturn in [Eu/Fe], vis¬ 
ible at low metallicities, is a consequence of the fixed amount of 
Eu produced by NSM, coupled with the paucity of NSM events 
and the constant mixing volume assumed in our model. When 
a NSM pollutes a simulated box early on, it produces a value in 
the [Eu/Fe] vs [Fe/H] space, dependent on the mass of the pre¬ 
vious enriching SNell. The volume enriched by NSM and SNII 
with the lowest amount of iron creates the upper tip of this up¬ 
turn towards low metallicity. Then in all the volumes polluted 
by NSM, the probability of having another Eu enrichment is 
low, so they evolve towards lower [Eu/Fe] and higher [Fe/H] 
by the subsequent enrichment of Fe by SNII, creating the diag¬ 
onal shape from high [Eu/Fe] with low [Fe/H] to low [Eu/Fe] 
with higher [Fe/H] shown in Fig [T] (left panel). A more detailed 


dynamical treatment, where the pollution in not confined to 
fixed volumes, does not produce this sharp feature dArgast et al.l 


2QQ4I) . Also to consider a variation in the r-process production 


(as in the models NS03-NS04, see below) leads to a smoothing 
of this sharp trend. The fraction of Eu-free stars, shown in the 
band at [Eu/Fe]=-2.4, is not negligible for this model and also 
for the next NSM models; we discuss the implication of this 
outcome in the Sect. 4.4 in more detail. 


A model with a delay time of 10 Myr (model NSOl) essen¬ 
tially behaves the same (see Fig.[T] central panel), while delay 
times as long as 100 Myr cannot be accepted, because they lead 
to even worse predictions, with all stars with Eu measurements 

below [Fe/H]-2 dex being unexplained by the model (model 

NS02, Fig.[T] right panel). Therefore, in the following we only 
discuss models that assume At^s= 1 Myr (apart for model 
NS-rMRD). Similar results have been found by Matteucci et 
al. (2014). 


In Fig. (left panel) we show the predictions of model 
NS03, which is the same as model NSOO, except that the 
masses ejected by NSMs are not constant, but instead ran¬ 
domly distributed around the mean value, 5-10“^ M© (see 
equation ID). This approach is similar to the one adopted by 
ICescutti & Chiappinif (l2014l) to explain the Ba, Y, Sr, and Eu 
abundance scatter in halo stars in the framework of a model 
where MRD SNe are the only r-process element sources. In 
Fig. [2] (central panel) we show the predictions of model NS04, 
which differs from Model NS03 in the adopted value of the 
average Eu yield from NSMs, 2.5-10“^ M©, i.e. half 

the value adopted in Model NS03. Moreover, Model NS04 as¬ 
sumes that the fraction of massive stars leading to NSMs is 
twice the val ue adopted in Model NS 03, namely 4%. It is worth 
noticing that iMatteucci et all (12014 ) require that 2% of mas¬ 
sive stars are hosted in binary systems, leading to NS merg¬ 
ers in order to fit the current merger rate of these systems in 
the Galaxy. However, there are no strong arguments against 
a variation in this quantity in the early Universe. As a matter 
of fact. Model NS04 leads to theoretical predictions that agree 
with the data much more closely than Model NS03. In partic¬ 
ular, Model NS03 also predicts a wide spread at intermediate 
metallicities, which also implies a non-negligible population of 
stars with sub-solar [Eu/Fe] ratios at [Fe/H]-2, which are ac¬ 

tually not observed. We note that, however. Model NS04 does 
not allow fully addressing the existence of stars with [Eu/Fe] <0 
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Fig. 1. Left panel, the results for [Eu/Fe] vs [Fe/H] for model NSOO. This model has a delay time for the neutron star mergers 
of 1 Myr, constant Eu production of 5-10“^ M© per merging event, and no Eu from SNell. The density plot is the distribution of 
simulated long-living stars for our model (see the bar below the figure for the colour scale). The long-living stars formed withou t 
Eu (formally [Eu/Fe]=-ex?) are shown at [Eu/Fe]=-2.4. The model predictions are compared to data collected in iFrebell ( 201Qh : 
we show as black dots stars with [C/Fe] <0.9 (to avoid binary enrichment), the open dots are stars with no carbon measurement. 
Central panel, same as left panel, but for model NSOl. This model has a delay time for the neutron star mergers of 10 Myr, 
constant Eu production of 5-10“^ M© per merging event, and no Eu from SNell. Right panel, again same as for the left panel, but 
for model NS02. This model has a delay time for the neutron star mergers of 100 Myr, constant Eu production of 5-10“^ M© per 
merging event, and no Eu from SNell. 


for [Fe/H] < -3. In Fig 12] (right panel), we show the results of 
Model NS03, but with a mass range for neutr on star progeni 
tors of 9-50 M©, instead of 9-30 M©. In fact, iMatteucci et al 


(I 2 OI 4 ) suggest that Eu in the Galaxy can only be reproduced 
by NSM if the range of progenitors of neutron stars extends up 
to 50 M©. However, this change does not provide significant 
variations in our predictions, mostly because of the stochas¬ 
tic nature of our model in which the enrichment at the lowest 
[Fe/H] is not necessarily due to the most massive stars. 

We are aware that to obtain a reasonable agreement 
with the observational data, we had to assume the shortest 
times cale for the neutron star s mergers suggested in the liter¬ 
ature ( Belczvnski et al.|[2002h . and this can be an extreme as¬ 
sumption. In fact, a progenitor lifetime of 1 Myr is already the 
sh ortest possible timescale in the distribution of merger times 
by iBelczvnski et all (l2002h : then, unfortunately, their popula¬ 
tion synthesis method is not very predictive of the merging 
timescale: “Unlike all other binary properties ... the qualita¬ 
tive characteristics of the merger-time distributions appear to be 
rather sensitive to a number of model parameters.” Moreover, 
the local rate of short-hard gamma ray bursts can provide at 
least a rough observational constraint on the typical lifetimes 
of its progenitors (so NS-NS binaries), and the results seem 


to indicate lifetimes of a few Gyrs (see iNakar et ^ l2006h . 
Finally, the estimated merger times for the fe w NS-NS bina¬ 
ries in our Galaxy are of the order of 100 Myr (iKalogera et al 


200 ll: IPiran & ShavivluOO^ . so all these observational indica¬ 
tions cannot exclude that some NSM explode on a very short 
timescale, but at least they seem to point to an average merger 
time that is much longer than 1 Myr. 


Argast et al.l (l2004l) have also studied the NSMs as a pos¬ 


sible source of r-process material. Their model dynamically 
traces the diffusion of the ejecta, so in this respect it is more 
detailed than our model. However, it is not clear how well their 
model follows the chemical evolution of the Galaxy, since no 
comparison with the global MDF of the halo was provided, and 
no outflow from the system was considered. Still, despite these 
differences, the results do not differ substantially, apart from the 
tendency of their model to produce too many r-process-poor 
stars at intermediate metallicity compared to our results and 
presumably a larger number of metal free stars. Also their con¬ 
clusions are similar: “NSM as major r-process sources are only 
consistent with observations” if “ the NSM population has coa¬ 
lescence timescales shorter than approximately 10 Myr.” They 
underline that NSM could be co-producers of the r-process ,to- 
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Fig. 2. Left panel, same as Fig.[Tl but for model NS03. This model has a delay time for the neutron star mergers of 1 Myr, variable 
Eu yield per merging event (average value 5-10“^ M©), and no Eu from SNell. Central panel, same as left panel, but for model NS04. 
This model has a delay time for the neutron star mergers of 1 Myr, variable Eu yield per merging event (average value 2.5-10“^ M©), 
and no Eu from SNell. In this model, the fraction of massive stars leading to NS merging events is twice what is assumed in previous 
models. Right panel, same as left panel, but we present here a model similar to model NS03 in which we consider a mass range for 
neutron stars progenitors of the NSM of 9-50 M© instead of 9-30 M©. 


gether with core-collapse SN, but they do not analyse the pos¬ 
sible outcome. We provide this comparison in the next section. 


4.2. NSM and core-collapse models 

A possible solution to the issue of the existence of stars with 
[Eu/Ee]<0 at [Ee/H]< -3 is to consider Eu production from 
core-collapse SNe. At the present time, a scenario in which all 
SNe II produce Eu is not supported by nucleosynthesis mod¬ 
els. A standard SNII explosion can produce chemical elements 
only up to the Sr-Y-Zr peak du ring a normal core-collapse SN 
explosion dArcones et al.l12007 ). Nevertheless, at least to have 
a possible idea of the impact of this production, we present 
the results of model NS-fSN in Eig.[3](left panel). This model 
shares the prescriptions for Eu production in NSMs with model 
NS03, but add to it a contribution from SNell increasing from 
2-10“^ M© for a 20 M© star to 5-10“^ M© for a 50 M© star. 
As for the synthesis of Eu, Model NS-fSN adopts prescrip- 
tions similar to those of Model Mod2SNNS in Matteucci et al.l 


(120141) . The model proved successful in reproducing the aver¬ 
age [Eu/Ee]-[Ee/H] relation of solar neighbourhood stars, the 
solar Eu abundance, and the present [Eu/H] gradient across the 
disc. 

Finally, we explore the possibility that a fraction of mas¬ 
sive stars end their life with a strong magnetic field and 
an extremely fast rotation in the ir inner cores. This M RP 
SN scenario ha s been studied in Iwinteler et al.l (2012 h and 
Nishimura et all ( 2015h . According to their results, the explo¬ 
sion is able to trigger a production of r-process material. We 


assume that this channel is active only in 10% of massive stars 
and only in the metal-poor environments (log(Z)< -3) since 
this channel is favoured at low metallicities. As in previous 
cases, in Model NS-i-MRD, the MRD production is coupled 
with the NSM productions, but with a long delay of 100 Myr. 
In Fig. [3] (right panel) the results of Model NS-hMRD are pre¬ 
sented. In this case the model is able to explain the distribution 
of the data at low metallicity and the overall trend very well. 

Model NS-hMRD provides a possible solution if the 
timescale for Eu enrichment from NSM turns out to be rel¬ 
atively long. Clearly, more detailed assumptions, such as a 
time delay distribution for the NSM or a fading of the MRD 
contribution (rather than a switch off at a certain metallicity), 
would be an improvement in the modelling. However, at the 
moment neither observational nor theoretical indications allow 
us to build better scenarios, so we prefer to keep things sim¬ 
ple. Nevertheless, we underline again that this double scenario 
is a promising way to explain the distribution of [Eu/Fe] in ex¬ 
tremely metal-poor stars. 

4.3. NSM and spinstar scenario 

Now we turn to other neutron-capture elements, with both an 
r- and an s-process component, namely Sr, Ba, Zr, and Rb. In 
the first column of Fig. [4] we show the resulting [Sr/Ba] vs 
[Fe/H] for Models NS-rMRD and NS03, compared to Model 
MRD-i-s B2. In Model MRD-i-s B2, we assumed that 10% of 
the massive stars end their lives as MRD SNe, and we con¬ 
sidered the possibility that the amount of mass ejected as r- 
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Fig. 3. Left panel: same as Fig.[Tl but fo r Model NS+SN . The prescriptions about Eu nucleosynthesis for Model NS+SN are 
similar to those of model Mod2SNNS bv iMatteucci et alJ ( 2014 ) , but it allows for inhomogeneous mixing in the early Galaxy. 
Right panel: same as left panel but for Model NS+MRD. In this model we consider the production of Eu by neutron star mergers 
with a delay of 100 Myr and also a production by MRD SN (10% of all SNe II) but only for Z < 10“^. 


proce ss varies following equation l4l(^see lCescutti & Chiappinil 
2014 for further details). The aim is to investigate the impact 
of the different r-process productio ns assumed in NS+M RD 
and NS03 on the previous r esults of ICescutti et alJ ( 2013 ) and 
Cescutti & Chiappinil ( 20141) . in particular on the idea of spin- 
star production as a viable explanation for the spread in [Sr/Ba] 
in the observat i onal d ata, which is a concept introduced by 
Chiappini et alJ (l201lh . We recall that in all the models, Sr, 
Ba, Zr, and Rb are produced by both s- and r-processes. The 
r-process yields are obtained by just scaling the Eu contri¬ 
bution according to th e chemical abundan ce ratios observed 
in r-process-rich stars dSneden et alJ l2008l) . We consider the 
contribution by spinstars , following the yields calculated by 
Erischknecht et al. ( 20121). T he original Model MRD+s B2 (see 
Cescutti & Chiappini 20141) was in fact able to reproduce the 
scatter of [Sr/Ba] data at [Ee/H]< -2.5, thanks to the combina¬ 
tion of the spinstar production with [Sr/Ba] >0 and the r-process 
production fixed at [Sr/Ba]—0.2. 

Model NS+MRD turns out to be very similar to Model 
MRD+s B2, and we note just tiny variations in the predicted 
distribution of long-living stars. This happens because the mod¬ 
els have essentially the same nucleosynthesis prescriptions at 
low metallicity (log(Z) < -3). At log(Z) > -3, the models 
are different: NSMs are the only r-process producers in Model 
NS+MRD, but the impact of this change is mitigated by an ISM 
that is already homogenous and quite rich in neutron capture 
elements, so the models do not show remarkable differences. 

On the other hand, an important difference arises in the 
comparison between the NS03 model and the other models. 
Since the NSM rate is lower than the assumed rate of MRD 
SN events, the average time needed to have an enrichment of 


the r-process by NSMs is longer in a stochastic realization. 
Therefore, we find in this model volumes that are only polluted 
by spinstars for a longer period, and they can survive without 
r-process pollution up to [Fe/H]—2. This can be noticed in the 
more extended prediction of long-living stars with high [Sr/Ba] 
and also in higher stellar density, in the area [Sr/Ba]>0 and 
-3 <[Fe/H]< -2. 

In the second and third columns of Fig. |4l we display our 
new results on the impact of spinstars in the chemical evolu¬ 
tion of [Zr/Ba] and [Rb/Ba]. We underline that these chemical 
predictions for Rb are the first results available for this ele¬ 
ment (probably due to the paucity of observational data). For 
Zr we have found few homogeneous inodels that describe it s 
evolution (among them iTravaglio et al]l2004 ICescuttil 12007 ). 
The model results for [Zr/Ba] are very similar to the results 
for [Sr/Ba]; since the observational data have similar distribu¬ 
tions, we can conclude that our models reproduce the stellar 
abundances well at low metallicity for Zr. Also, the differences 
amongst the predictions of the different models, in particular 
between Model NS03 and Models MRD+s B2 and NS+MRD, 
follow the same pattern for the [Sr/Ba] case. In contrast, the 
predictions for [Rb/Ba] are more complex and show different 
trends. The spinstars for these elements do not always produce 
a ratio above the r-process signature, but in both directions 
(above and below). 

The different behaviour of Rb with respect to Sr (and Zr) 
has several explanations. The first is that the ratio [Rb/Ba] is 
generally lower than [Sr/Ba]. This is because Rb is less stable 
than Sr (and Ba). This first difference is enhanced in spinstars 
for the following reasons. The isotopes produced in spinstars 
that are most abundant have different properties for Sr and Rb, 
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Fig. 4. From the left [Sr/Ba], [Zr/Ba] and [Rb/Ba] vs [Fe/H] in the halo; the density plot i s the distribut ion of simulated lon g-living 
stars fo r our models. Superimposed, we show the abundance ratios for halo stars (data from iFrebe J20 1 (ih . The symbols for the lFrebeJ 
(1201 Ol) data are black dots for normal stars, and black open circles for stars without carbon measurement. 


in particular their positio n along the s-process path and their 
neutron-richness (see e. g. lKappeler et al ■1201 ll and references 
therein). For Sr the most abundant isotope produced is ^^Sr, 
which is the most neutron-rich stable isotope of Sr. ^^Sr is a 
bottleneck of the s-process path at N = 50, and it is one of 
the main isotopes produced by the weak s process in massive 
stars. Its production is generally high. The situation is differ¬ 
ent for Rb. This element has only two stable isotopes, ^^Rb 
and ^^Rb (as opposed to 4 for Sr). The Rb isotope most abun¬ 


dantly produced by the weak s process is ^^Rb. The production 
of ^^Rb is sensitive to the branching point at ^^Kr. If neutron 
densities are high, more ^^Kr capture another neutron instead 
of beta-decaying to ^^Rb. This thus reduces the production of 
^^Rb significantly. This reduction is compensated for partially 
by an increase in the production of ^^Rb (via beta decay of 
^^Kr), but generally the production of ^^Rb is an order of mag¬ 
nitude lower than ^^Rb. Neutron densities are high in the spin- 
star models considered in this work (with a low y) rate). 
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This leads to a generally lower production of Rb relative to Sr 
for the reasons listed above and explains why [Rb/Ba] are lower 
than [Sr/Ba] (and [Zr/Ba]). 

Indeed, the different behaviour of Rb with respect to Sr (and 
Zr) is an interesting signature that can be investigated in fu¬ 
ture observational campaigns; at the moment, only a few mea¬ 
surements at higher metallicity in globular clusters are avail¬ 
able. Again, the NS03 model predicts a stronger presence of 
long-living stars polluted by spinstars compared to the models 
MRD-rs B2 and NS-fMRD. 

In summary, the assumptions on the r-process production 
of the models NS-fMR D a nd NS03 do not alter t he fin dings 
of ICescutti et all (12013h and ICescutti & Chiappinil (l2014l) . and 
the spinstars produce a spread in [Sr/Ba] and [Zr/Ba] for these 
models matching the observational data. For [Rb/Ba], no data 
are available at the moment to check our results, but we pre¬ 
dict a different spread than for [Sr/Ba] and [Zr/Ba]. The NS03 
model tends to enhance the population of long-living stars pol¬ 
luted by spinstars alone, compared to the two other models. 
This point is not easily tested, but the tendency does not appear 
to be di splayed at least in the availa ble observational data (see 
Fig. 6 in ICescutti & Chiappinil2Q14 ). since stars with very high 
[Sr/Ba] (> 0.5 dex) are not so frequent. 


4.4. Eu- free stars? 



[Fe/H] 


Fig. 5. Ratio of Eu-free stars over the total number of 
stars for bins of 0.5 dex in [Fe/H]. The model results 
are displayed with a black dashed line (model NS04), red 
solid line (model NS-i-MRD), and blue dotted line (model 
NS-f SN). The red squares are the the observational proxy 
for this ratio, so the ratio between stars in which Eu only 
presents an upper limit over the number of stars for which 
at least Ba has been measured. The error bars are only 
plotted to show the dimension of each bin in [Fe/H]. 


The different r-process sites analysed in the previous sec¬ 
tions (MRD SN, NSM) have a common feature: only a small 
fraction (<10%) of the massive stars actually enter Eu pro¬ 
duction. On one hand, this promotes the observed spread of 


[Eu/Fe] at low metallicities. On the other hand, it ensures that 
at extremely low metallicities a complementary fraction of low 
mass stars are formed in regions where no r-process synthesis 
took place; for this reason, it seems likely to form some Eu-free 
stars. The possibility that the r-process production is delayed 
further increases the fraction of Eu-free stars, since more stars 
can be formed before any Eu enriches the ISM. 

Our stochastic models predict the formation of a fraction 
of Eu-free stars: all the low mass stars formed in each stochas¬ 
tic realization before the first enrichment by a r-process site 
are Eu-free. Therefore, we investigate how this result matches 
the observations of Galactic halo stars. We check this by com¬ 
puting the ratio of Eu-free stars over the total number of stars 
in our models in a given [Ee/H] bin and then compare these 
theoretical ratios to an observational proxy for this ratio. The 
number of observed stars where it has been measured is only 
an upper limit for Eu compared to the total number of stars for 
which Ba has been observed (see Eig. B- We are aware that 
there is probably a bias in this plot, since it is affected by the 
weakness of the Eu line, so below a certain metallicity, the S/N 
to distinguish the Eu line cannot be achieved during standard 
observations. Eor this reason, we expect this proxy to only be 
an upper limit of the real fraction of Eu free stars. 

Model NS-i-SN always shows a lower fraction of Eu-free 
stars compared to our observational proxy; since given the 
above, this proxy is an upper limit, it is a positive result. The 
model NS -I- MRD predicts a higher number of Eu-free stars than 
in our observational proxy. Nevertheless, a slightly higher rate 
for the MRD events at extremely low metallicity could solve 
this problem. Model NS04, which is the one performing better 
in the pure NSM scenario, always predicts a too high ratio of 
Eu-free stars compared to our observational proxy, and it also 
only drops to zero at very high metallicity. 


We note that also rotating massive stars can produce a tiny 
amount of Eu by s-process (^see lErischknecht et alJ2ni2h . how- 
ever the calculations at the present time can produce at best 

an [Eu/Ee]-4 at [Ee/H]-3. We have indeed included the 

yields for spinstars in the models discussed in this paper, but 
their impact is simply too low to be seen in the plots. 


Should future observations show that Eu is always present 
in the composition of EMP stars, similar to Ba and Sr (Roederer 
2013), the most likely explanation would be that ah core- 
collapse SNe are producing at least a tiny fraction of Eu 
through some r-process channel (as in our model NS-l-SN). 
Another way to solve the riddle, as investigated by Komiya et 
al. (2014), would be to assume that the observed Eu is accreted 
on the stellar surfa ce from the ISM, but this result contradicts 
Erebel et ^ (l2009h 


Overall, the results obtained with the inhomogeneous 
model for the chemical evolution of the Galactic halo adopted 
and reinforce previous conclusions by 
) that either NS mergers explode on a 
very short timescale or that at least a fraction of SNell (MRD 
SNe) must produce some Eu. At present, the last scenario 
should be preferred, since it is the one that best reproduces the 
available observations. 


m this work conhrm 


Matteucci et all (2014 
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5. Conclusions 


References 


In this paper, we adopt the inhomogeneous mo del for the chem 


ical evolution of the Galac tic halo presented bvICescuttil ( 2008 ) 


i 


2014 


and further developed in ICescutti & Chiappinil (201 
and ICescutti et alJ (12013b . to study the evolution of the n 
capture elements Eu, Sr, Ba, Zr, and Rb in the early Galaxy. 
We implemented the Eu production from co alescing neutron 


stars in the code following the prescriptions of iMatteucci et al 
( 2014 ). These in turn rest on recent work bv iKorobkin et al 


(12012h showing that large amounts of r-process elements, 0.01 
Mo, can be produced by NS when they merge. Out of these, the 
mas s of newly produced Eu can lie in the range -10“^ - 10“^ 
Mq. iMatteucci et all (12014 have already studied the effect of 
such a large Eu production from NSMs by means of a detailed 
homogeneous chemical evolution model for the Milky Way. 
They suggest that, though NSMs can in principle explain the 
history of Eu enrichment in the Milky Way by themselves, a 
more realistic situation requires that both NSMs and (at least) 
a fraction of core-collapse SNe contri bute to the Eu produc- 
tion in the Milky Way. However, while IMatteucci et al.l (12014 
could explain the observed trend of [Eu/Ee] versus [Ee/H] in 
the solar vicinity, the abundance of Eu measured in the Sun 
and the present-day abundance gradient of Eu along the disc, 
they did not address the problem of explaining the wide spread 
of [Eu/Ee] ratios measured in halo stars, because their analysis 
was carried out with the use of homogeneous chemical evolu¬ 
tion models, which are intended to reproduce average trends. 

In the present work we provide a complementary analysis 
by using stochastic chemical evolution models and investigate 
the information contained in the scatter of several abundance 
ratios (e. g. [Eu/Ee], [Sr/Ba], [Zr/Ba], and [Rb/Ba]). We show 
that either NSMs explode with a very short fixed timescale or 
both channels - NSMs and core-collapse SNe (at least a frac¬ 
tion of them as in the MRD scenario) - must be active. In order 
to explain the spread, as well as both the presence of stars with 
low [Eu/Ee] at [Ee/H] < -3 and a certain fraction of Eu free 
stars, the last scenario should be preferred. 

The resu lts shown here are con sisten t with the recent con¬ 
clusio ns by ICescutti et al.l ( 2013 ) and ICescutti & Chiappinil 
( 2Q14 of an important role of spinstars in the early Universe, 
not only as contributors to light elements, but also to n-capture 
elements such as Sr, Ba, Zr, and Rb. Indeed, their inclusion, 
together with NSM and MRD SN sites, leads to models that 
can reproduce the scatter observed in all of these elements very 
well. 
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